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Radial stability of the continuous pressure gravastar is studied using the conventional 
Chandrasekhar method. The equation of state for the static gravastar solutions is derived and 
Einstein equations for small perturbations around the equilibrium are solved as an eigenvalue 
problem for radial pulsations. Within the model there exist a set of parameters leading to a 

cN : 

stable fundamental mode, thus proving radial stability of the continuous pressure gravastar. 
It is also shown that the central energy density possesses an extremum in p c (R) curve which 



a* 



represents a splitting point between stable and unstable gravastar configurations. As such 



in 

the p c (R) curve for the gravastar mimics the famous M(R) curve for a polytrope. Together 
with the former axial stability calculations this work completes the stability problem of the 
continuous pressure gravastar. 
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I. INTRODUCTION 

m 
> 

(f-) ■ Gravitational collapse as the stellar nuclear fuel is consumed could lead to black holes - objects 

in , 

£T) , which are accepted by scientific community but their undesired and even paradoxical features 

(singularities, horizon) have motivated research in a direction of finding massive objects (stars) 
without singularities and without horizon. One of these alternatives is a gravastar. 

Since the seminal work of Mazur and Mottola [1;] the concept of the gravitational vacuum star - 
the gravastar - as an alternative to a black hole has attracted a plethora of interest. In this version 
of the gravastar a multilayered structure has been introduced: from the repulsive de Sitter core 
(where a negative pressure helps balance the collapsing matter) one crosses multiple layers (shells) 
and without encountering a horizon one eventually reaches the (pressureless) exterior Schwarzschild 

pr n 

spacetime. Later some simplifications [2|, 13J and modifications [4| have been introduced in the 
original (multi)layer - onion-like picture. 

An important step was done when it was shown that due to anisotropy of matter comprising the 
gravastar [5(] one can eliminate layer(s) and, by a continuous stress-energy tensor, the transition 
from the interior de Sitter spacetime segment to the exterior Schwarzschild spacetime is possible 
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fl] (see also [J). 



(see also [7(). The pressure anisotropy in the spherically symmetric geometry was perhaps 
first introduced by G. Lemaitre g|] and suggested by Einstein (as quoted in [8]). The vanishing 
radial pressure with transversal pressure only was shown to be enough to support a stable object. 
Further development has brought different refinements to the original anisotropy notion. 

The pressure anisotropy, which is shown to be a necessary condition for the existence of a gravastar 



,5, 



is met also in boson star models 



111 ] and wormholes 12] . The anisotropy (defined as a difference 



between the transversal and radial pressure) vanishes at the center (r = 0) of the star as well as at 

nn 

its boundaryfr = it). The gravastar has been confronted with its rivals - black holes [la. 1141 ] and 



wormholes 12 . 



and investigated with respect to energy conditions (violations) [17J] and its 
3, Ql- An interesting question has been posed several times: is it possible to 



charged properties 
distinguish the gravastar from a black hole 



20M22I]? In Ref. [20] it was shown that gravitational 



radiation could be used to tell a gravastar from a black hole. However, the definite answer to this 
question has not been given at satisfactory level and the gravastar research is still a dynamic field 



with recent papers like fllgJl etc. 

Almost every research mentioned above to some extent addresses the problem of the gravastar 
stability, since the stability problem is crucial for any object or situation to be considered as 
physically viable. In Ref. 1] it was first shown that such an object is thermodynamically stable, 
while axial stability of thin-shells gravastars was tested in j], 4\. Stability within the thin shell 
approach based on the Darmois- Israel formalism was recently reviewed in [^j- I n 2Q] stability 



analysis of the thin shell gravastar problem is closely related to an attempt to distinguish the 
gravastar from a black hole by analysis of quasi-normal modes produced by jaxial perturbations. 
Problem of stability of a rotating the thin shell gravastar was addressed in 
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3l| 



(multi)layer version of the gravastar was also considered in [14I . H . 

The axial stability of the continuous pressure gravastar was shown to be valid in 
analysis was based on the Ref. 



2$. Stability in the 



This 



321 ] where stability of objects with de Sitter centers was investigated. 



In this paper we analyze the radial stability of the continuous pressure gravastars [6] following 



the conventional Chandrasekhar method. Origina 
testing the radial stability of the isotropic spheres 



ly Chandrasekhar developed the method for 



331 ] in terms of the radial pulsations. In Ref. 



3a ] Chandrasekhar's method was generalized to anisotropic spheres. Stability of anisotropic stars 



was investigated before in 



37 



and radial stability analysis for anisotropic stars using the quasi- 



local equation of state was given in 39[]. The standard mathematical procedure is here applied to 
an object with a peculiar behavior of pressures (see below) and although the mathematical rigor 
was never abandoned, the analysis due to the character of the object could be considered as a toy 
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model analysis of radial stability. 

This paper is organized as follows. In the next section the linearization of the Einstein equations 
is given. Static solutions are described, an equation of state is derived and the pulsation equation 
is obtained. In Sec. Ill the eigenvalue problem for the radial displacement is presented. Results 
and discussion are given in the last section. 

Unless stated explicitly we shall work in units where Gn = 1 = c. 

II. LINEARIZATION OF THE EINSTEIN EQUATION 

In this paper the response of the continuous pressure gravastar model to small radial pertur- 
bations is considered. Assuming that the pulsating object retains its spherical symmetry, one can 
introduce the Schwarzschild coordinates: 

ds 2 = e^dt 2 - e x ^dr 2 - r 2 d9 2 - r 2 sin 2 6d<p 2 , (1) 

where A and v are, in this dynamical setting, time-dependent metric functions. 

The standard anisotropic energy-momentum tensor appropriate to describe continuous pressure 

gr avast ars is: 

T; = (p + Pr )u^ - g»p r + y (p* - Pr ) + k^ipt ~ Pr), (2) 

where vP is the fluid 4- velocity, vP = dx^/ds, In and kn are the unit 4- vectors in the 9 and 4> 

directions, respectively, = —r 5^, = —rsmOSf,. 

The velocity of the fluid element in the radial direction £ is defined by: 

where £ is the radial displacement of the fluid element, r — > r + £(r, t). The components of the 
4- velocity are obtained by employing u^u^ = 1 and Eq. ([3]): 

u x = ( e -^ 2 ,e e - y / 2 ,0,0). (4) 

The non-zero components of the energy-momentum tensor ([2]) linear in £ are: 

Tt = P , T r r =-p r , T$ = T* = -p u T t r = i(p + p r ), Tt = e x -»£(p+p r ). (5) 
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The components of the Einstein tensor for the metric JT]) are: 

1 
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e- ( X - - 


\ r 
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1A 1 (7) 

(8) 
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— + - + — + + -e"^ I A + — - — . (9) 
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Following the standard Chandrasekhar method, all matter and metric functions should only slightly 
deviate from their equilibrium solutions, 

A(r, t) = A (r) + 5X(r, t), u(r, t) = u (r) + Su(r, t), (10) 

p(r,t) = p (r) + 5p(r,t), p r (r,t) = p r o(r) + Sp r (r,t),p t (r,t) = p t o(r) +5p t {r,t). (11) 

The subscript denotes the equilibrium functions and Sf(r,t) are the so-called Eulerian pertur- 
bations, where / G {A, v, p,p r ,pt, }. The Eulerian perturbations measure a local departure from 
equilibrium in contrast to the Lagrangian perturbations, denoted as df(r,t), which measure a 
departure from equilibrium in the co- moving system (fluid rest frame). The Lagrangian perturba- 
tions in the linear approximation play a role of a total differential and are linked to the Eulerian 
perturbations via the equation (see e.g. Ref. [401]): 

df(r,t) = 6f(r,t) + ti(r)t;. (12) 

A linearization of the Einstein equations G 1 ^ = 8ttT^ leads to the two sets of equations: one for the 
equilibrium (static) functions and the other for the perturbed functions. The equilibrium functions 
obey the following set of equations: 

+ ^ ( 13 ) 

4 + 7 + T + 4 (15) 

In practice, one usually combines these three equations into the Tolman-Oppenheimer-Volkoff 
(TOV) equation: 

1 2 

PrO = (PO +Pr0) "a + ~^0, (16) 
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where Ilo denotes the anisotropic term IIo = pto—PrO- The other set of equations emerging from the 
linearization of the above Einstein equations yield the set of equations for the perturbed functions: 



- x "5X"j =8irr 2 5p, (17) 
bv = (i/ + ^\ SX + 8Trre Xo 5p r , (18) 



e -Ao 



SX = -87r£Go +Pro), (19) 



^°- Uo (po + Pro)i + \{po + ProW + \{b~P + 5p r )v' Q + 5p' r - hu = 0. (20) 



361 ] and it serves to probe the radial stability 



Equation (|20p is known as the pulsation equation 
of the system of interest. It is actually the TOV equation for the perturbed functions which is 
obtained - analogously as the non-perturbed TOV - by combining Einstein equations for perturbed 
functions. 

In order to solve the pulsation equation (|20p for the gravastar all perturbed functions should be 
expressed in terms of the radial displacement £ (and its derivatives) and the equilibrium functions. 



In performing this, one first integrates Eq. (119p in time, yielding: 

SX = -87rre Xo C(po+Pro). (21) 
Using this expression for 5p in Eq. (I17p one obtains: 

Sp = -^[r 2 (po+Pro)^}' ■ (22) 

After inserting 5X in Eq. (|18p a dependence on 5p r remains, which should be expressed in terms of 
the displacement function (and its derivatives) and the equilibrium functions. To accomplish this, 
one ought to explore the system at hand in more detail. 



33] 



One of the possibilities, as suggested firstly by Chandrasekhar for isotropic structures 
and more recently by Dev and Gleiser for anisotropic objects 361 ] . is to make use of the baryon 



density conservation to express the radial pressure perturbation in terms of the displacement func- 
tion and the static solutions. In this approach the adiabatic index appears as a free parameter. 
Chandrasekhar used this method to establish limiting values of the adiabatic index leading to an 
(un)stable isotropic object of a constant energy density. He showed that there were no stable 
stars of this kind if the adiabatic index was less than 4/3 + kM/R (k is a constant of order unity 
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depending on the structure of the star, M and R are the star's mass and radius). In Ref. [36] 
the Chandrasekhar method was extended to various anisotropic star models and showed that the 
limiting value of the adiabatic index is shifted to lower values, i.e. anisotropic stars can approach 
the stability region with smaller adiabatic index than in the Chandrasekhar's case. 

In this paper our primary concern is to probe the radial stability of one particular anisotropic 
object - the gravastar. Due to the peculiar character of the gravastar (especially its radial pressure 
- see below) one cannot expect the adiabatic index to be constant along the whole object. In fact 
the adiabatic index is a function of the energy density and pressure(s). This is the main reason why 
in this paper stability will not be tested by fixing the appropriate values of the adiabatic index that 
guarantee stability. The required information will rather be extracted from a given static solution 
by constructing the equation of state. 

A. Static solution 



The procedure discussed so far is applicable to all spherically symmetric structures. To ap- 
ply it to gravastar configurations one has to recall the basic characteristics of gravastars in the 
continuous pressure picture Q]. The energy density po(r) is positive and monotonically decreases 
from the center to the surface. Gravastars have a de Sitter-like interior, p r o(0) = —po(0), and a 
Schwarzschild-like exterior. Furthermore, the atmosphere of the gravastar is defined as an outer 
region, near to the surface, where "normal" physics is valid [fj, i.e. where both the energy-density 
and the radial pressure are positive and monotonically decreasing functions of the radius. In the 
gravastar 's atmosphere the sound velocity v s , with 



2 _ d PrO 



(23) 



dp ' 

is real (v^ > 0) and subluminal (v s < 1). 

From the peculiar shape of the gravastar's (radial) pressure one can immediately infer that the 
sound velocity ought to be real only in the gravastar's atmosphere, whilst in the gravastar's interior 
it is imaginary, < 0. This is the main reason why, in probing the radial stability, we shall be 
primarily concerned with the physical processes occurring in the gravastar's atmosphere. 

To construct a static gravastar, the energy density profile and the anisotropic term are adopted 
from the previous work [a, : 

Po(r) = p c (l-{r/R) n ), (24) 
n (r)=/3p (0"Vo(r). (25) 



7 




FIG. 1: The energy density po/pc, radial pressure p r o/Pc tangential pressure pto/Pc and compactness po 
as a function of radius r/R for {R,n,m} ~ {1,2,3}. Three different values of the central energy density 
p c = {0.190,0.202,0.210} and their anisotropy strengths /3 = {92.905,81.410,76.110} correspond to the 
lower, middle and upper curve, respectively, ro denotes the radius at which the sound velocity (|23|) vanishes 
(for the central curve). 

Here n, m are (free) parameters and p c = po(0) is the central energy density. /3 is the 
anisotropy strength measure and R is the radius of the gravastar for which p r Q(R) = 0. po(r) 
is the compactness function defined by Ho(r) = 2mo(r)/r, where mo(r) is the mass function 
mo(r) = Air J po(r)r 2 dr. The radial pressure p T o is a solution of the TOV (fl~6|) and the tan- 
gential pressure is readily obtained from the anisotropy and the radial pressure by employing the 
identity pto = p r o + IIo- The particular form of the anisotropic term is dictated by the behavior of 
pressures at de Sitter core, since at r = the anisotropy should vanish as seen from (16). Also, 
the above form of the anisotropy term ensures that the radial pressure vanishes at r = R. The 
transversal pressure vanishes as well although it is not necessary to be the case (see Ref. [6J for 
the gravastar model with non vanishing transverse pressure.). The anisotropy strength measure is 
controlled as well by the energy conditions which have to be met. 

One such solution for fixed (R,n,m) = (1,2,3) is shown in Fig. 1 for three different values 
of the central energy density p c corresponding to three different values of the anisotropy strength 
(3. Since the radius R is fixed there is an interplay between the central energy density p c and 
anisotropy strength f3 - a higher central energy density p c requires a smaller anisotropy strength /3. 
We shall elaborate on this particular choice of parameters in Section IV, where the radial stability 
of these three gravastar configurations will be tested. 

In the inset of Fig. 1 the radial pressure close to the surface is extracted in order to show important 
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features of the gravastar's atmosphere. At the radius ro the sound velocity of the fluid vanishes 
(dpro / dpo\r=r = 0) and hence ro serves as a division point of the propagating (or physically 
reasonable) (r > ro, Vg > 0) and non-propagating regions (r < ro, v 2 s < 0) when probing radial 
pulsations of the gravastar. 

The dominant energy condition (DEC), i.e. p r o,Pto < po, is obeyed by both radial and tangential 
pressure throughout the gravastar. The compactness function has also been shown in Fig. 1. 

B. Equation of state 

In this subsection we note that the equation of state (EoS) appropriate to describe the 
gravastar (inferred from the input functions (|24p and (|25p) is actually a functional of the energy 
density (only), parameterized by the anisotropy strength f3. Next this result is used to compute 
the Eulerian perturbation of the radial pressure 5p r from the EoS, by perturbing the energy 
density only. Ultimately this completes the task to express all perturbed functions in terms of the 
displacement (and its derivatives) and the static solutions. 

Generally, for isotropic structures, before solving the TOV, one assumes that the pressure p and 
the energy density p are functions of the specific entropy s and the baryon density n. If a system 
is described by the one-fluid model, then in static and dynamic settings it exhibits isentropic 
behavior (constant s), in which case one can set s = 0. Thus it is possible to eliminate the baryon 
density n and express the pressure in terms of the energy density only, leading to a barotropic 
equation of state p = p(p). It is a rather simple task now to perturb this EoS and express the 
perturbed pressure in terms of the perturbed energy density. 

For anisotropic objects the EoS is highly dependent on the anisotropic term model (see e.g. 
the TOV p6p ). The particular choice of the anisotropic term used here (j25jl is a functional (or 
a quasi-local variable 1 ) of the energy density. This means that for a fixed anisotropy strength 
j3 there is a two parameter family of values {p c ,R} belonging to the same EoS (see Fig. 2). As 
a consequence, one can obtain perturbed (radial) pressure by perturbing the energy density only, 
and keeping the anisotropy strength (3 fixed. 

To illustrate this in more detail let us introduce an analytic form of the EoS which, to a good 

1 By the quasi-local variable we mean a function which is an integral in space of some local function - for example, the 
mass function mo(r) is a quasi-local variable of the energy density (which is a local function) as it is the volume 
integral of the energy density (the same holds for the compactness function). For a discussion of quasi-local 
variables and quasi-local EoS see e.g. Refs. |3^] and Ref. [3S| | . 
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FIG. 2: The radial pressure p r o/Pc is plotted against the energy density po/pc (EoS) for {R, n, m} — {1,2,3}. 
Three different values of the central energy density p c = {0.190, 0.202, 0.210} and their anisotropy strengths 
j3 = {92.905, 81.410, 76.110} correspond to the lower, middle and upper curve, respectively. 

approximation, describes the gravastar configuration defined by (i24"|) and (l25j) : 2 

ProlPo) = ~Po (y~ a tolPol) ■ ( 26 ) 

Here a is closely related to the anisotropy strength (3, p>o[po] is the compactness function which is 
a functional of the energy density. Now it is clear that for a fixed a the (radial) pressure is fully 
determined by the energy density. 

Hence, following the reasoning outlined above and making use of Eq. (|12p . in the linear approx- 
imation the Eulerian perturbation for the radial pressure is: 

6 P r = -p' r o( + ^j^-(6p + P , oO- (27) 



Here dp ro [po]/dpo denotes functional derivative of the radial pressure with respect to the energy 

dpro/dr 
dpo/dr 



density. This is equal to d j ro /j T as both the radial pressure and the energy density are functions 



of radius r only. 

Similarly, the Eulerian perturbation of the anisotropy 5U assumes the form: 

5U = -U' C + ^(6p + p'S. (28) 
apo 



With the above two expressions the pulsation equation (|20j) is fully determined. However, before 
proceeding to solve the pulsation equation it is useful to rewrite Eq. (|27p in a slightly different form 



2 It is worth noting that the analytic form of the EoS (|26p is not restricted to the chosen energy density (|24[) . For 

_ 2 

example, it is also appropriate to describe a gravastar with the energy density of the form po{r) = p c e 17 r . 
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in order to compare the result in this paper with that of Chandrasekhar's for isotropic, and Dev 
and Gleiser's for anisotropic stars. By means of the TOV (I16p the perturbed energy density (I22D 
can be written as 

5P = -Pot ~ (Po+Pro)— (rV^) - -H £. (29) 
Inserting this result in Eq. (|27|) the radial pressure perturbation becomes 

S Pr = -PU ~ (Po + Pro) d -^^ (rV-^Y - 2n %M e . (30) 
If one now identifies the adiabatic "index" as 

P0 + PrO dpro[po] , . 

7 = - A , (31) 

PrO dpo 



361 ] . Eq. (86). However, the 



the result derived in Eq. (|30j) reduces to that of Dev and Gleiser 
expressions for 7 differ. Moreover, if one turns off anisotropy (Ho = 0) Chandrasekhar's result is 
obtained. 



III. THE PULSATION EQUATION AS AN EIGENVALUE PROBLEM 

As in the Chandrasekhar method all matter and metric functions exhibit oscillatory behavior 
in time, f(r,t) = e lujt f(r). Hence the pulsation equation assumes the form: 

V e + Vx? + V 2 t = -uj 2 VU, (32) 

where VoiViiVi and are polynomial functions of r, depending on the static solutions only (see 
Fig. 3) . Eq. (|32p represents an eigenvalue equation for the radial displacement £ (with ui 2 being an 
eigenvalue). Solutions of this differential equation are obtained by specifying boundary conditions 
in the center and at the surface of the gravastar: 

£ = at r = 0, (33) 



Ap r 



at r = R. 



(34) 



The boundary condition in the center demands that there is no displacement of the fluid in the 
center of the gravastar. The boundary condition at the surface follows from the requirement 
that the Lagrangian radial pressure perturbation has to vanish at the surface 40H42] . In the 



model presented here where Ap r = (dp r o[po]/dpo)Ap, the sound velocity vanishes at the surface, 



11 




FIG. 3: The polynomial functions To, Vi/20, P2/IOO and Vu from the pulsation equation (f2"U|) for {R, n, m} = 
{1,2,3} and {p c ,P} = {0.202,81.410}. 

dp r o/dpo\ r= R = 0. This means that, apart from being finite, there are no further restrictions on 
Ap(R). This also implies that it is sufficient to demand that £(R) and £'(R) are bounded in order 



to satisfy the boundary condition at the surface 



4l|. The choice 



?(R) = 0, (35) 

enables one to compare the results in the gravastar's atmosphere with the radial oscillations of the 
polytropes. This can be relevant as the EoS of the gravastar's atmosphere close to the surface can 
be approximated by the polytropic EoS p r oc p 1+1 ' n p, where n p is a polytropic index 
In order to study radial stability of the system described by Eq. ([20]) subject to the boundary 



conditions (f33|) and flMJ), it is plausible to recast the pulsation equation into the standard Sturm- 
Liouville form (see e.g. Ref. [401]): 

(vey + = -^ 2 w£, (36) 

where 

v = e JPi/v dr and Q = ^lr, w = ^J±v. (37) 

Po Po 

The leading coefficient in the pulsation equation Vq has three zeros - two at the ends {0, R} and 
one in the interior region ro (dp r o/dpo\ r=ro = 0), hence Vi/Vq has three singular points (see Fig. 



3), though all three are regular singular points or Fuchsian singularities 



m 3 - 



3 A singular point r* is regular (or Fuchsian) if the function Vx/Vo has a pole of at most first order, and the function 
V2/V0 has a pole of at most second order at the singular point r = r*. 
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In order to obtain V the integral J Vi/Vo dr should be calculated, and since the interior singu- 
larity arises at tq which is a division point between propagating and non-propagating domains, it 
is reasonable to divide the whole interval Z = (0, R) in two parts: X\ = (0, ro) and I2 = (ro,R). 
In performing the integration numerically infinitesimally small regions around all three singular 
points {0, ro,i?} are excluded, so that both integrals are rendered convergent and finite. As a 
consequence, the leading coefficient in the Sturm-Liouville equation V is a positive function on 
the (whole) interval X, whilst the weight function W is negative on the interval T\ and positive 
on the interval 12- As elucidated in the previous section, the interesting region is the gravastar's 
atmosphere, i.e. the second interval, X2. In this region the standard Sturm-Liouville eigenvalue 
problem formalism (see e.g. [411 ]) is applied, since V > and W > 0. Therefore if uj 2 is positive, 
co itself is real and the solution is oscillatory. If on the other hand uj 2 is negative, uj is imaginary 
and the solution is exponentially growing or decaying in time, thus signalizing instabilities. The 
number of nodes of the eigenvector £ for a given eigenvalue uj 2 is closely related to the stability 
criteria. To be more precise, if for uj 2 = eigenvector £ has no nodes, then all higher frequency 
radial modes are stable. Otherwise, if for uj 2 = eigenvector £ exhibits nodes, then all radial modes 
are unstable. Furthermore, if the system is stable, then the following relations hold 

ujI < w? < • • • < uj 2 n < . . . , (38) 

where n equals the number of nodes. 



IV. RESULTS AND DISCUSSION 



In testing stability of certain configurations in general, it seems natural that one attempts 
to find critical values of the parameters for which the system is marginally stable. Marginal 
stability means here that there exists a set of parameters for which the system exhibits the stable 
fundamental mode (n = 0) for Uq = 0. Then for the given set of parameters all higher frequency 
modes are radially stable. For example, in the case of neutron stars (described by the polytropic 
EoS), there exists a critical value of the central energy density for which the stellar mass M as a 
function of radius R is extremal. For such a critical value of the central energy density the star 
exhibits stable fundamental mode with ujq = 0, which implies that all higher frequency modes 
with the given central energy density are then radially stable. Furthermore, at the account of the 
M(R) curve one can then read off which configuration of the EoS will produce a stable star and 
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FIG. 4: The central energy density p c plotted against the radius R. For {R, n, m} — {1, 2, 3} the anisotropy- 
parameter f3 — {92.905,81.410,76.110} is constant on each curve and fixed by choosing the central energy 
density to be p c = {0.190,0.202,0.210} from the lower to the upper curve, respectively. The minimum of 
each curve represents marginally stable configurations. 

which will not. 

The continuous pressure gravastar model described here displays a quite similar behavior. For 
each EoS (fixed j3) the extremum of the p c {R) curve represents critical values of the parameters 
{p c ,R} for which the system exhibits a stable fundamental mode, Uq = (see Fig. 4). Then 
for such critical set of parameters all higher frequency modes are radially stable. Moreover, for 
smaller radii the system exhibits stability, whereas for larger radii (than the critical one) it reveals 
instability (see Fig. 4). In this sense the p c (R) curve for the gravastar mimics the well known M(R) 
curve for a polytrope. To prove these statements, the behavior of the displacement function £ is 
shown in Fig. 5 for three different cases. The radius R is, for simplicity, fixed (for all three cases) 
to be the critical radius of the central curve in Fig. 4. According to Fig. 4 one then expects that 
the radial displacements £ derived from the lower, central and upper curve in Fig. 4 will generate 
stable, marginally stable and unstable fundamental mode respectively. This is exactly what is 
shown in Fig. 5. The central (solid) curve in Fig. 5 represents the marginally stable fundamental 
mode - an eigenvector £ is obtained for Wq = 0. The upper (short-dashed) curve clearly shows 
stability of all radial modes as for Wq = there are no nodes, while the lower (long-dashed) curve 
reveals instabilities of all radial modes as there is a node in the fundamental mode. The lower, 
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FIG. 5: The displacement function £(r) for {R,n,m} = {1,2,3} and lu 2 = 0. Three different val- 
ues of the central energy density p c = {0.190,0.202,0.210} and their respective anisotropy strengths 
j3 = {92.905,81.410,76.110} correspond to the lower (unstable), middle (marginally stable) and upper 
(stable) curve respectively. 

middle and upper curves in Fig. 5 correspond to the upper, middle and lower curves in Fig. 1 
and Fig. 2, respectively. Here again one is able to relate this result to that of Ref. 36|]: from 



Fig. 5, according to the values of the anisotropy strengths j3, one can conclude that the anisotropy 
enhances stability. 

Albeit from the viewpoint of radial pulsations, the gravastar's inner region does not seem to be 
physically attractive as the sound velocity is imaginary there, it is important to add a couple of 
comments on the radial displacement's behavior in that region. Pulsations are strongly attenuated 
in the gravastar's interior (see Fig. 5). This holds for all uj 1 > 0. Therefore the radial pulsations of 
the gravastar as a whole can be seen as occurring prevalently in the gravastar's atmosphere whereas 
entering the interior region they are exponentially (but smoothly) attenuated. This is actually what 
one would intuitively expect from the repulsive gravitation caused by the de Sitter-like interior. 4 
In this paper the focus was set on one very specific star model - the gravastar. Therefore 
standard stability analysis which has been applied here in every detail could be considered as a toy 
model of the radial stability analysis. However it could be extended to a broader class of anisotropic 



4 A good example of such a space is an inflationary universe. The electric and magnetic fields of free photons in such 
an inflationary (quasi-de Sitter) space get (exponentially) damped as oc 1/a 2 , while the physical wavelength gets 
stretched as oc a. Here a denotes the scale factor of the Universe, which during inflation grows nearly exponentially 
in time. 
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stars with the anisotropy being a functional of the energy density. In this way the adiabatic "index" 
does not have to be set to a constant but calculated from the static configurations. This comprises 
one of the main results of this paper. 

The main result of this work is the observation that the continuous pressure gravastar-model 
presented here exhibits radial stability as illustrated in Fig. 4 and Fig. 5. This result is important as 
it, along with the axial stability analysis, suggests that gravastars, although not yet fully understood 
at the fundamental level, may be viable physical compact object candidates. 
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